Quantum Reversibility: Is there an Echo? 
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We study the possibility to undo the quantum mechanical evolution in a time reversal experi- 
ment. The naive expectation, as reflected in the common terminology ("Loschmidt echo"), is that 
maximum compensation results if the reversed dynamics extends to the same time as the forward 
evolution. We challenge this belief, and demonstrate that the time tr for maximum return probabil- 
ity is in general shorter. We find that tr depends on A = Eevoi/eprep, being the ratio of the error in 
setting the parameters (fields) for the time reversed evolution to the perturbation which is involved 
in the preparation process. Our results should be observable in spin-echo experiments where the 
dynamical irreversibility of quantum phases is measured. 



In this Letter we study the probabihty of return 
<2) for a generalized wavepacket dynamics scenario. 
The system is prepared in some initial state , which 
can be regarded as the outcome of a preparation pro- 
cedure which is governed by a Hamiltonian Tipjop- We 
assume that the quantum mechanical evolution is gener- 
ated by Hamiltonians with classically chaotic limit: The 
state is propagated for a time ti using a Hamiltonian 
Til ) and then the evolution is time- reversed for a time t2 
using a perturbed Hamiltonian 7^2 • The corresponding 
evolution operators are Ui and U2- The probability of 
return to the initial state is 



P{h,t2) = |(^p„p|C/2(t2)-^C/l(il)|*p.ep)|' 



(1) 



There are two special cases that have been extensively 
studied in the literature. The traditional wavepacket dy- 
namics scenario is obtained if we set t2 = 0. In this 
context the "survival probability" is defined as 



Psnit) = P{t,0) 



(2) 



The "Loschmidt echo" (LE) scenario is obtained if we set 
ti = t2 = t. In this context the "fidelity" is defined as 



(3) 



The theory of the fidelity was the subject of intensive 
studies during the last 3 years |illlilElllli,|i[l3- 

It has been adopted as a standard measure for quantum 
reversibility following and its study was further moti- 
vated by the realization that it is related to the analysis 
of dephasing in mesoscopic systems \vi\ . 

In the present Letter we consider the full scenario of 
a time reversal experiment. The probability to find the 
system in its original state is P{t) = P{t, 0) before the 
time reversal {t < T/2), and 



P{t) = P{T/2,t-T/2) 



(4) 



after the time reversal (< > T/2). The period T is the 
total time of the experiment. The naive expectation, 
which is also reflected in the term "Loschmidt echo", is 



to have a maximum for P(t) at the time t = T. We are 
going to show that this expectation is wrong. We find 
that the maximum return probability is obtained at a 
time tr which in general is shorter than that. Namely, 



T/2 < U < T. 



(5) 



If we have tr — T/2 we say that there is no reversibility. 
If we have tr T we say that we have a nearly perfect 
echo. We show that t^/T is a function of a dimensionless 
parameter < A < oo. Namely, 



I = /(A) 



A = ecvoi/ep 



(6) 



where £p,.ep quantifies the difference between the evolution 
Hamiltonian Ti. and the preparation Hamiltonian 7ip,,,p, 
while £„^oi quantifies the difference between the two in- 
stances Til and 7^2 of the evolution Hamiltonian, which 
are used for the forward and for the time-reversed evolu- 
tion respectively. The idea is that there is no way to have 
a complete control over the parameters (fields) of the sys- 
tems. Therefore there is an unavoidable difference (eovoi) 
between these two instances of Ti, which by the setup 
of the experiment are regarded as identical. The scal- 
ing function © takes the limiting value /(O) — 1 (echo) 
while for A > A* we get /(A) = 0.5 (no reversibihty). 
Here A* is some system-specific constant of order unity. 

The most popular preparation which is considered in 
the literature, either in the context of wavepacket dynam- 
ics or fidelity (LE) studies, is a Gaussian wavepacket. 
Obviously the choice of such preparation is motivated 
mainly by the wishful thinking of theoreticians. However, 
in many applications, one is not so much interested in 
evolving an initial Gaussian wavepacket. This is certainly 
the situation in quantum information processing [l^ and 
in spin-echo experiments lOj where one starts with a ran- 
dom initial state. Formally a Gaussian wavepacket can 
be regarded as the ground state of a phase-space shifted 
Harmonic oscillator. Therefore it is characterized by a 
very large £prcp, leading to A ^ 1. In this Letter we do 
not assume A ^ 1, but rather consider the general case. 
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In order to develop a general theory we need a model 
in which we have control over both ep,„p and £ovoi- 
We consider a quantized system whose classical ana- 
log has positive Lyapunov exponent. Its Hamiltonian 
H, = TC{Q, P;x) depends on a parameter (field) vari- 
able X which is determined by the experimental setup. 
For example it can be either a gate voltage or a mag- 
netic flux. The dynamics takes place within a classically 
small (but quantum mechanically large) energy window. 
The classical dynamics is assumed to have a well defined 
finite correlation time Td- We consider classically small 
(but possibly quantum mechanically large) perturbations 
(x I— > X -f Sx) . Accordingly, the Hamiltonian can be lin- 
earized as follows: 

H = £ + SxB (7) 

We define Hi and H2 by setting Sx = ±e„^„i. The re- 
quirement of having classically small Sx means that the 
phase space structure of Tii and Ti.2 is similar, and that 
any (small) difference in the chaoticity can be neglected, 
(we have verified that this smallness condition is satisfied 
for the example below). 

The preparation issue requires further discussion. 
The traditional possibility is to prepare a Gaussian 
wavepacket (also known as a coherent state preparation). 
We regard this possibility as uncontrolled because the 
value of Epjop is ill defined. To have a physically meaning- 
ful definition of ep,ep the natural procedure is as follows: 
We define a preparation Hamiltonian TYp^ep by setting 
Sx = Eprep. Then we start with an initial eigenfunction 
of £, and evolve it with ?ip„p until we get an ergodic-like 
steady state within an energy shell (note ^3|)- The width 
of this energy shell is proportional to ep„p. The resulting 
wavepacket is used as an initial state for the time reversal 
experiment. For the purpose of comparison we shall con- 
sider also a Gaussian wavepacket preparation. For such 
preparation A ^ 1 irrespective of its energy width. The 
reason is that this preparation does not occupy ergodi- 
cally its energy shell. Formally one can say that for a 
Gaussian wavepacket Hp^p differs enormously from the 
evolution Hamiltonian. There is no point in quantifying 
this difference. This is the reason why we use our con- 
trolled preparation procedure, where 7Yp„p differs from 
the evolution Hamiltonian in a well defined manner. In 
any case we shall verify that a Gaussian wavepacket is 
indeed like taking a preparation with A <g; 1. 

In our numerical investigation we use the model Hamil- 
tonian 

H{Q,P;x) = \{Pl+Pi + Ql+Ql) + xQlQl (8) 

with X = 1 4- Sx. It describes the motion of a particle 
in a 2D well (2DW). The physical units are chosen so as 
to have dimensionless variables. Therefore upon quan- 
tization the Planck constant ?i is a dimensionless quan- 
tity. Our numerical study is focused on an energy window 
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FIG. 1: Left panels: The probability P(t) to find the sys- 
tem in its original state, for the 2DW model. The two pan- 
els, from top to bottom, are for the representative values 
Ecvoi = 0.00315 and 0.26. We use various preparations: a 
coherent state preparation (solid line), a random superposi- 
tion with the same energy width (o), a random superposition 
(dot-dashed line) with a larger A where tpjop = 0.4, and a 
A > 1 preparation (dashed line) which was obtained by set- 
ting Epjcp = 0. The arrows indicate the time tr for the case 
Eprop = 0.4. Right panel: The simulations are now done with 
the corresponding RMT model. 



around « 3 where the motion is mainly chaotic with 
characteristic correlation time t^i « 1 13]. The quan- 
tization is done with h — 0.012. We write the Hamil- 
tonian matrix as in Eq. Q, using a basis such that £ 
is diagonal. The mean level spacing is A « 4.3 x Ti^ . 
As expected, on the basis of a general "quantum chaos" 
argumentation '':\4\, the matrix S is a banded matrix. 
More details regarding the band profile can be found in 
Ref. 0] . The only additional piece of information that is 
needed for the following analysis is the parametric scale 
Sxc- This is defined as the Sx which is needed in order 
to mix neighboring levels. It is given by the ratio A/cr, 
where a is the root mean square value of the near di- 
agonal matrix elements of the B matrix. For the above 
model Sxc « 3.8 * 

Fig. 1 (left panel) displays representative results of 
simulated time reversal experiments. One experiment is 
done with a coherent state preparation, and we indeed see 
behavior that looks like an echo (t^ ^ T). Qualitatively 
the same behavior is observed for a random preparation 
that has A <C 1 . Once we take a preparation with a larger 
value of A, we realize that the compensation time is in 
general tr < T. In particular with the A > 1 preparation 
we do not observe any quantum reversibility (tr — T/2). 

In Fig. 2 we present some results for tr/T. We clearly 
see that tr is smaller for larger values of Eevoi • Much more 
illuminating is Fig. 3 where we present our numerical 
results for tr and various T. The points corresponding to 
the same A but different values of e^^^i and 6^,^^ fall onto 
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the same smooth curve with a good accuracy, confirming 
the scahng hypothesis ©. Moreover, we see that for 
A > A*, with A* « 1.3, we have tr = T/2 (no echo) 
irrespective of the actual value of ep„p and Eevoi- Thus, if 
we want a A > A* type of preparation, we simply can set 
£prop — 0, which means an eigenstate of £. From Fig. 3 
we also see that in the other limiting case of A ^ 1 we 
have (nearly) an echo. This is the same as in the case of 
a Gaussian wavepacket preparation. 

In order to explain the observed results for tr we look 
on Fig. 4. There are two axes along which the theory 
of P{ti,t2) is quite well known. One is the traditional 
wavepacket dynamics axis t2 = along which Psn{t) is 
defined, while the other is the "LE axis" (^2 = ti) along 
which P^Eit) is defined. In the figure we also indicate the 
course of a time reversal experiment. It is clear that in 
order to have tr < T the contor line P(ii, ^2) = -Psr(2^/2) 
should meet the ti axis in a sharp angle. Furthermore, 
if we want to have some remnant of an echo at the end 
of the period {t — T) we have to cross the contour line 
P{ti,t2) = Ple(T/2). The condition for that is 

Psr(T/2) < P^e{T/2) (9) 

In the following discussion we would like to assume that 
both £prop and Cevoi are larger than JcCc, which means that 
the perturbations are strong enough to mix levels. In 
such case the decay of Psnit) or Pi^sit) is approximately 
exponential: 

Psn{t), P^^it) = exp(-7i) (10) 
Moreover in both cases 7 is given by the expression 

7 = min(7pT,7sc) (H) 



where 7pt is the value which is determined by pertur- 
bation theory, while 7sc is determined by semiclassical 
considerations. 

Once details are concerned the theory behind the ex- 
ponential approximation Eq. H1(J|I becomes quite different 
in the two respective cases [PsR{t), Pi^E{t)]- The theory 
of the survival probability is related to the parametric 
theory of the LDOS [l3lll5| . Namely, for relatively small 
perturbations 7 = 7pt oc (ep^p/^Xc)^ is essentially the 
width of Wigner's Lorentzian, while for large perturba- 
tions 7 = 7sc oc Epjep is the width of the energy shell (in 
the latter case the exponential approximation is at best 
a good fit). In contrast to that, the theory of the fidelity 
Pi^sit) is related to a theory of dynamical correlations, 
and cannot be reduced to the LDOS analysis 0. The 
best theory to date is semiclassical (see |9| and refer- 
ences therein) . As in the case of the survival probability 
we have 7 = 7pt oc (s^^ai/Sxc)'^ for small perturbations, 
while for large perturbations 7 = 7sc, in some typical 
cases, is related to the Lyapunov exponent. We would 
like to point out that the study of the general condi- 
tions for having a fingerprint of the Lyapunov exponent 
in time reversal experiments is still an open issue for fu- 
ture study . The existing semiclassical theory for the 
"echo" phenomena assumes Gaussian wavepackets. 

By definition Psr{T/2) depends only on £p„p, while 
Ple{T/2) is mainly sensitive to Eovoi- Therefore it is evi- 
dent that a small A is a condition for having < T. If we 
have A < 1 then the contour line P(ii,t2) = Pi^e{T/2) 
can be very close to the LE-axis. This implies that for 
A <C 1 we can get nearly a perfect echo behavior {tr ^ T) . 
This picture, as we have seen before, is supported by our 
numerical findings. 
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FIG. 2: The time tr is calculated from the simulation of P(t), 
and the ratio tr/T is presented for various values of the period 
T, and for various values (see legend) of £ovoi. An average over 



a number of initial wavepackets with the same £p 



0.4 was 



performed. We clearly see that tr is smaller for larger values 
of Eovoi. The left panel corresponds to the 2DW simulations, 
while the right panel to the effective RMT model. 




FIG. 3: The ratio tr/T for different initial random wavepack- 
ets which axe determined by Eprop, and for various Ecvoi. The 
horizontal axis is the scaled variable A. The data scale nicely 
in accordance with Eq. JSJ. The dashed line is a simple poly- 
nomial fit. The left panel corresponds to the 2DW simulations 
while the right one to the corresponding RMT model. 
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FIG. 4: The contour line of P{t\,t2) that goes through the 
point (r/2,r/2). One solid line ("no echo") is for the case 
A > 1. The other solid line is for the case of a relatively 
small A. The dashed line illustrates the course of a time re- 
versal experiment, while the "LE axis" is the line along which 
Ple is defined. 

As we have seen above, the applicability of semiclassi- 
cal considerations is not essential for having an "echo". 
The general picture that we have outlined should be valid 
also in the absence of a semiclassical limit. This is in 
contrast to the impression that one might get from the 
recent literature. In order to establish this provocative 
statement, in a way that leaves no doubts, we use a simple 
random matrix theory (RMT) procedure. We take the 
resulting banded matrix B of the 2DW model ||HJ), and 
randomize the signs of the off-diagonal terms. In this 
way we get an effective RMT (ERMT) model of the typ e 
that had been introduced by Wigner 50 years ago |l7| . 
The model is characterized by the same mean level spac- 
ing, and by the same band-profile as the physical 2DW 
model. Consequently the generated dynamics is charac- 
terized by the same correlation time (the latter is deter- 
mined by the bandwidth). But unlike the 2DW model, 
the ERMT model is lacking a semiclassical limit. In the 
right panels of figures 1-3 we demonstrate the results of 
simulations that were done with the ERMT model. We 
clearly see that we get similar results (in Fig. 3 the RMT 
drop is slightly sharper). 

In summary, we have shed a new light on the physics 
of quantum reversibility, and in particular we have in- 
troduced the concept of compensation time tr, which re- 
places the misleading terminology of "echo" . Our predic- 
tions should be tested in wave field evolution experiments 
such as s pin polarization echoes in nuclear magnetic res- 
onances [lOl [iSj. In particular we have considered the 
realistic case of a general preparation, and clarified the 
role of semiclassical considerations in the theory. 
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